Probing modifications of General Relativity using current cosmological observations 
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We test General Relativity (GR) using current cosmological data: the cosmic microwave back- 
ground (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs- Wolfe (ISW) effect from 
the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation 
of supernovae Type la (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the 
weak lensing (Wb) data from CFHThS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear 
and mildly non-linear scales. We first test a model where the effective Newton's constant, n, and 
the ratio of the two gravitational potentials, 77, transit from the GR value to another constant at 
late times; in this case, we find that GR is fully consistent with the combined data. The strongest 
constraint comes from the ISW effect which would arise from this gravitational transition; the ob- 
served ISW signal imposes a tight constraint on a combination of /i and r\ that characterizes the 
lensing potential. Next, we consider four pixels in time and space for each function /1 and 77, and 
perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes 
are consistent with GR within the errors. Only one eigenmode shows a 2<r deviation from the GR 
prediction, which is likely to be due to a systematic effect. However, the detection of such a devia- 
tion demonstrates the power of our time- and scale-dependent PCA methodology when combining 
observations of structure formation and expansion history to test GR. 



I. INTRODUCTION 

As cosmological observations improve, new possibili- 
ties arise for testing the physics that governs the evo- 
lution of our Universe. Precise all-sky measurements of 
the CMB by the WMAP satellite [1] have established 
that cosmic structure developed from a nearly scale- 
invariant initial spectrum of adiabatic fluctuations [2]. 
Baryon acoustic oscillations from SDSS [3] and grow- 
ing catalogs of supernovae [4], in combination with the 
CMB, have tightened the constraints on the background 
expansion history, indicating a strong preference for the 
cosmological concordance model, ACDM. Correlating the 
CMB anisotropies from WMAP with wide-sky catalogs 
of galaxy counts have made it possible to detect the ISW 
effect [5, 6], obtaining independent evidence for the ac- 
celerating expansion of the Universe [7, 8] . Weak Lensing 
measurements by surveys such as COSMOS [9, 10] and 
CFHTLS [11, 12] have afforded the use of shear correla- 
tion functions and power spectra in order to test cosmol- 
ogy. 

In parallel, the problem of cosmic acceleration has 
motivated explorations of new theoretical ideas, includ- 
ing the possibility that GR may be modified on large 
scales. Anticipating the substantial improvement in cos- 
mological data sets that is expected with surveys such as 
DES [13], Pan-STARRS [14], LSST [15] and Euclid [16], 
model-independent frameworks for testing GR against 
observations of the growth of cosmic structures have been 
developed [17-26]. Recently, there has been progress to- 
wards a consensus on what properties such a framework 



should have [21, 25, 26]. First, it should be general, i.e. it 
should be able to describe a wide range of modified grav- 
ity models. Second, it should not violate the consistency 
of super-horizon sized perturbations with the background 
expansion. Third, it should involve as few parameters as 
possible, while still being flexible enough to capture most 
of the significant information contained in the data. 

A good pragmatic starting point is to look for evidence 
of departures in the fundamental relationships among the 
perturbative fields familiar in cosmology: matter density 
and velocity perturbations, and the metric perturbations 
[19, 27-29]. It is commonly agreed that one needs two 
general functions of scale and time, in addition to the 
conservation of energy-momentum, to specify the evolu- 
tion of the linear perturbations. For example, one can 
introduce fi, relating the gravitational potential to the 
density contrast, and 77, which relates the gravitational 
potential to the spatial curvature. Given these functions, 
it is possible to calculate the cosmological perturbations 
and all the observables with a full Boltzmann integration 
code, such as MGCAMB [25, 26]. 

The actual parametrization of the functions has not yet 
been researched as fully. The goal is to strike a balance 
between simplicity, i.e. working with as few parameters as 
possible, and allowing for enough flexibility in these func- 
tions to capture all of the information contained in the 
data. For example, one can discretize /j(fc, z) and rj(k, z) 
on a grid in (fc, z) space and treat their values at each 
grid point, which we will call pixels, as independent pa- 
rameters [30]. Motivated by simplicity, one might wish 
to start with the simplest possible model, such as con- 
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stant values of /i and 77, and add complexity only if the 
fitted parameter values of this simple model show hints 
of departure from ACDM. However, this logic may not 
work in modified gravity studies and might lead to miss- 
ing important information contained in the data. As the 
PCA analysis in [30] shows, the shapes of the well con- 
strained eigenmodes suggest a higher sensitivity to scale- 
dependent features in /z and 77, compared to their average 
values or the time dependence. 

In this work, we first use a Fisher forecast-based PCA 
to determine the minimum number of pixels necessary 
to describe the shape of the best constrained eigenmodes 
of /1 and 77. We then use a Monte Carlo Markov Chain 
(MCMC) algorithm to fit these parameters to a combina- 
tion of the available data, including CMB, ISW, SNe and 
WL, and find constraints on their de-correlated combina- 
tions. Throughout the paper we assume that the back- 
ground expansion is given by the flat ACDM model which 
is strongly favored by the current constraints on the ex- 
pansion history, and look for deviations from its predic- 
tions for the density perturbations. From a theoretical 
perspective, flatness is motivated by the inflationary ori- 
gin of the Universe, and the viable models of modified 
gravity studied in the literature tend to be indistinguish- 
able from ACDM at the background level. 

Other tests of GR have been performed, in which dif- 
ferent choices were made for the functions /x and 77 (or 
a related set of parameters). In [31, 32] they were taken 
to have a specific form of time dependence, while in [33] 
they have a specific form of time and scale dependence; 
finally, [32] allowed them to vary in three redshift bins. 
The results of these studies show a good consistency with 
ACDM. To compare with the results of some of these 
studies, we consider a single-transition-in-redshift model 
in addition to our scale-dependent PCA method, gener- 
ally finding a good agreement with GR. In addition to al- 
lowing for scale-dependence, other important differences 
between our study and the treatment in [32] include us- 
ing the ISW cross-correlation data, using the CFHTLS 
WL data coming only from linear and mildly non-linear 
scales, and simultaneously varying two functions /x and 77, 
while [32] varied them one at a time when working with 
a 3-bin model. The key conceptual difference from [33] 
is that we do not use WL data from a deeply non-linear 
regime, and we do not assume a specific scale and time 
dependence of the functions /x and 77, but rather perform 
a PCA of their values on a gird in (k, z). 

We find that the agreement with ACDM is statisti- 
cally more significant when no scale-dependence is al- 
lowed. After performing a PCA, we find that seven of 
the eight eigenmodes are consistent with GR within the 
errors. One eigenmode shows a 2a deviation from the GR 
prediction, but can be directly traced to a feature in the 
WL aperture mass dispersion spectrum at 120 arc min, 
which is most likely caused by a systematic [34]. However, 
the detection of this effect shows the benefits of adopting 
a more flexible scale-dependent pixcllation of /x and 77, 
and demonstrates that using scale-independent methods 



could potentially hinder the detection of new physics, or, 
as in this case, simply the better understanding of the 
data. 

The paper is organized as follows. Section II reviews 
the data sets used in this work. After describing our 
parametrization in Section III, we present the constraints 
on departures from GR in Section IV, and finally draw 
conclusions in Section V. 



II. OBSERVABLES AND DATA 

In this section, we summarize the observables that will 
be used to constrain deviations from GR, and explain the 
data sets for these observables. 



A. Integrated Sachs- Wolfe effect 

The ISW effect [5] is a secondary anisotropy of the 
CMB which is created whenever the gravitational po- 
tentials are evolving in time. This is due to the net en- 
ergy gain that the CMB photons acquire when traveling 
through varying potential wells, and it is therefore a di- 
rect probe of the derivatives of the potentials $, "P. In 
more detail, this effect generates additional CMB tem- 
perature anisotropies in any direction n given by 

9 IS w(n) = ^P^(n) = + [r, n(r - r)] dr, 

J- isw J K ' 

(1) 

where t is the conformal time, the dot represents a con- 
formal time derivative and the integral is calculated along 
the line of sight of the photon. 

A direct measurement of this effect is difficult, due to 
the overlap with the primary CMB anisotropies, whose 
amplitude is at least 10 times bigger. An additional prob- 
lem is that the ISW signal is biggest on the largest an- 
gular scales, which are most affected by cosmic variance. 
It is nevertheless possible to detect this signal by cross- 
correlating the full CMB with some tracers of the large- 
scale structure (LSS) of the Universe [6]: the primary 
CMB signal, generated at early times, is expected to have 
null correlation with the LSS, while the ISW anisotropies, 
produced at low redshift, correlate with the LSS distribu- 
tions since they trace the fluctuations in the potentials. 

We can then use a galaxy survey with visibility func- 
tion dN/dz(z) as a tracer of the LSS, and we can write 
the galaxy density fluctuation in a direction rii as 



<*s(ni) = 




where b g is the galaxy bias and S rn the matter density 
perturbation. On the other hand, the ISW temperature 
anisotropy in a direction n 2 is given by 

eisw(n 2 ) = - J A ($ + *) (n 2 , z)dz, (3) 
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where e _K ( z ) is the photons' visibility function and n the 
optical depth. After choosing a particular data set for 
the CMB and the LSS, we can then define the auto- and 
cross-correlation functions as 



(0) = (0(^)5^)) 
'(0) = (6^)5^)), 



(4) 
(5) 



where is the full CMB temperature anisotropy and 
the averages are calculated over all pairs at an angular 
separation $ = |ni — n 2 1 • Alternatively, the above calcu- 
lation can be written in harmonic space, and the auto- 
and cross-power spectra are then derived. 

We use the ISW data from [7] , which were obtained by 
cross-correlating multiple galaxy catalogs with the CMB 
maps from WMAP. The data used trace the distribution 
of the LSS in various bands of the electromagnetic spec- 
trum, with median redshifts 0.1 < z < 1.5, and consist of 
six catalogs (infrared 2MASS, visible SDSS main galax- 
ies, luminous red galaxies and quasars, radio NVSS, and 
X-ray HEAO). This is an approximation of a true tomo- 
graphic study of the ISW signal. 

All maps were pixellated on the sphere, with a pixel 
size of 0.9 deg. The measurements were done in real 
space, calculating the angular cross-correlation functions 
(CCFs) between the maps. These were linearly binned in 
steps of 1 deg for angles deg < $ < 12 deg; so the data 
set consists of 78 points (cj 9 ) h s - 

A well known property of the correlation functions is 
that their data points arc highly correlated; in this case, 
in particular, the high degree of correlation is present also 
between data points belonging to different catalogs, due 
to the partial overlaps in redshift and in sky coverage of 
the sources. For this reason, the full covariance matrix 
between all data points C, 3 - is a very important piece 
of information, and it was estimated in [7] using several 
Monte Carlo and jack-knife methods. Here we use the 
matrix produced with the most complete technique, a 
full Monte Carlo method where both galaxies and CMB 
maps were simulated and then correlated to measure the 
expected noise and covariance. 

The calculation of the likelihood of a particular model 
given the ISW data is done as follows. First, the the- 
oretical CCFs (cJ' ff )theo and auto-correlation functions 
(ACFs) are calculated with a full Boltzmann integra- 
tion within MGCAMB, based on the redshift distributions 
of the sources. The galaxy bias parameters are assumed 
to be independent and constant for each catalog, and 
are rescaled for each model imposing that the ACFs on 
small angular scales match the observations. Finally, the 
theoretical CCFs are multiplied by this rescaled bias to 
calculate the Xisw distribution, given by 



Xisw 



J2 [(cJ 9 Us - (cf 9 )thco] [C- 1 ] 



x 



B. Supernovae and Cosmic Microwave Background 

For the SNe data, we use the sample combination la- 
beled (e) shown in table 4 of [4], which is a compilation 
of the SDSS-II SNe sample plus Nearby SNe, ESSENCE, 
SNLS and HST. To calculate the SNe likelihood, we use 
values from the MLCS2K2 light curve fitter, and marginal- 
ize over the nuisance parameter, which is the calibration 
uncertainty in measuring the supernova intrinsic magni- 
tude. Note that [4] found a discrepancy of the constraints 
on the FwCDM model (standard CDM model in a flat 
Universe plus a dark energy component with a constant 
w) using the two SNe fitters MLCS2K2 and SALT-II. How- 
ever, the discrepancy is much smaller (within la) for a 
flat ACDM model as shown in Tables 13 and 17 in [4], 
namely, 

fi m = 0.312 ± 0.022(stat) ± 0.001 (syst) (MLCS2K2) 
O m = 0.279 ± 0.019(stat) ± 0.017(syst) (SALT - II) 

Since we will assume the background evolution is the 
same as that in the flat ACDM model, the choice of the 
SNe fitter does not affect our final results significantly, 
but we should bear in mind that systematic errors are 
now comparable to statistical errors in SNe observations. 

For the Cosmic Microwave Background data in our 
analysis, we use the WMAP five-year data including the 
temperature and polarization power spectra [38, 39], and 
calculate the likelihood using the routine supplied by the 
WMAP team 1 . 



C. Weak lensing 

We use the cosmic shear observations from the 
CFHTLS-Widc third year data release T0003 [11, 35], 
in which about 2 x 10 6 galaxies with ZAB-magnitudes be- 
tween 21.5 and 24.5 were imaged on 57 sq. deg. (35 sq. 
deg. effective area). We use the aperture- mass disper- 
sion M ap [36] following [11, 12]. As in those studies, one 
can obtain the relevant % 2 by fitting the theory-predicted 
aperture- mass dispersion (M^tiico given by a model pa- 
rameter vector p to (M^ p ) b s measured at angular scales 



XMap 



(P) = E(( M ap(^)) obs -<Ma 2 p(^,P)) thoo )[C- 1 ] y 
i 

(«(%))ob S -< M a 2 p(^,P)) thco ). (7) 



'J 
X 



Note that the data covariance matrix C is the one used 
in [11] and [12]; it contains shape noise, (non-Gaussian) 
cosmic variance and residual B-modes [34] . Since it is dif- 
ficult to model the weak lensing nonlinearity in modified 



(cj 9 )thoo 



(0) 



1 http://lambda.gsfc.nasa.gov/ 
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gravity in a model- independent way [37] , we only use the 
aperture-mass dispersion data measured between 30 and 
230 arc min, to remove the strongly nonlinear region from 
the data. For angles smaller than 30', the difference be- 
tween linear and non-linear predictions becomes greater 
than a factor of two, and we do not wish to suppose that 
non-linear corrections are reliable on smaller scales. 

Theoretically, the aperture-mass dispersion is related 
to the weak lensing power spectrum via [36] 



(M 2 ap )(6) 



dee 

~2T 



P K (£) 



24J 4 (6>£) 



(8) 



for the choice of filter in [11]; here, the lensing power 
spectrum P K is a projection of the 3D matter-density 
power spectrum P$, weighted by the source galaxy red- 
shift distribution and geometric factors, and J Q (x) is the 
Bessel function of the first kind. To model the redshift 
distribution of the galaxies, we follow [11] and use the 
parametrization 



n{z) cc 



z a + z ab 



Z» + C 



n(z)dz=l, (9) 



where J\f D {a, b, c} is a set of nuisance parameters to be 
marginalized over, and we have imposed Gaussian priors 
on them following [11]: a = 0.612 ± 0.043, b = 8.125 ± 
0.871, c = 0.620 ± 0.065. The distribution is normalized 
by setting z max = 6. Then we can calculate the x 2 f° r 
the redshift uncertainty as, 



E 



m - n\Zi 



(10) 



;th 



where rii is the normalized number of galaxies in the 
redshift bin and n(zi) the fitting function, evaluated at 
the center of the redshift bin. As described in [11], the 
uncertainty <7j of n, contains Poisson noise, photo- z error 
and cosmic variance, and we neglect the cross-correlation 
between different bins. Then we obtain the \ 2 f° r weak 
lensing in the same way as [11], 



XWL 



2 i 2 



(11) 



D. Further priors 

Finally, we impose la Gaussian priors on the Hubble 
parameter and baryon density of h — 0.742 ± 0.036 and 
Ob/i 2 = 0.022 ± 0.002 from the measurements of Hubble 
Space Telescope (HST) [40] and Big Bang Nucleosynthe- 
sis [41] respectively, and a top hat prior on the cosmic age 
of 10 Gyr < t < 20 Gyr. The total likelihood is taken 
to be the product of the separate likelihoods C of each 
dataset we used; thus the total \ 2 is the sum of separate 
X 2 from individual observations plus that from the priors 
if we define x 2 = ~ 2 log£. 



III. THE PARAMETRIZATION OF MODIFIED 
GROWTH 

In order to test gravity against the growth of cosmo- 
logical perturbations in a model-independent way, one 
needs a generalized set of equations to evolve linear per- 
turbations without assuming GR. We work within the 
framework of the Boltzmann integrator MGCAMB [25] 2 , 
which is a variant of CAMB [45] 3 . This is based on a sys- 
tem of equations that allows for a general modification of 
gravity at linear order in the perturbations, while respect- 
ing the consistency of the dynamics of long-wavelength 
perturbations with the background expansion [43, 44]. 
An interested reader can find a detailed discussion of 
the equations used in MGCAMB and their comparison with 
other methods in the literature in [26]. 

We consider scalar metric perturbations about a FRW 
background for which the line element in the conformal 
Newtonian gauge reads 



ds 2 



\t) [(1 + 2*) dr 2 - (1 - 2$) dx 2 



(12) 



where $ and VP are functions of time and space. We as- 
sume adiabatic initial conditions and covariant conser- 
vation of the energy-momentum tensor of matter. The 
matter conservation at linear order provides two equa- 
tions which in Fourier space can be written as 



5' + ^-v-W = 0, 
an 

k 

v' + v -<£ = , 

an 



(13) 
(14) 



where 5 is the energy density contrast, v the irrotational 
component of the peculiar velocity, and primes indicate 
derivatives with respect to In a. In order to solve for the 
evolution of the four scalar perturbations {8, v, we 
need two additional equations, provided by a theory of 
gravity (such as GR), which specify how the metric per- 
turbations relate to each other, and how they are sourced 
by the perturbations of the energy-momentum tensor. 
One can parametrize these relations as 



— = r](a,k), 



k 2 ^ 



-47rGa 2 fi(a,k)pA , 



(15) 
(16) 



where A is the gauge- invariant comoving density contrast 
defined as 



A = 6 



aH 



(17) 



r?(a, k) = n(a, k) = 1 in GR, while in an alternative model 
/i and rj can in general be functions of both time and 
scale [24, 25, 42]. 



2 http : //userweb . port . ac . uk/ - zhaog/MGCAMB . html 
1 http://camb.info/ 
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Defining fi and 77 in this particular way makes 
Eqs. (13)-(16) consistent on all linear scales. As shown 
in [26], on super- horizon scales \i naturally becomes ir- 
relevant and we are left with only one function, 77, as 
expected from the super-horizon consistency conditions 
[43, 44]. Also, having /j defined through the Poisson equa- 
tion involving A, as opposed to 5, allows for fi to be equal 
to unity on all scales for GR. 

From {77, fi}, we can derive other parameters which 
may be more suitable for interpreting observational con- 
straints. For example, since we are using WL and ISW 
observations in this paper, we will be essentially measur- 
ing the power spectra of the lensing potential ($ + \E r ) and 
its time derivative. On the other hand, 77 is not probed 
directly by any observable and would be highly degener- 
ate with /1, as also pointed out by [32]. For that reason, 
in addition to {77, /i}, we also present our results in terms 
of another function, E, defined as 



E(a, k) 



fc 2 (* + $) _ /_t(l + rj) 
8nGpa 2 A ~ 2 



(18) 



Note that specifying fj, and £ is equivalent to working 
with fj, and 77. We use both parametrizations to discuss 
the physics and to interpret our final results. 

Since we arc interested in testing GR at late times, we 
assume /_i(a, k) = rj(a, k) = 1 at early times, namely, for 
z > z s where z s denotes the threshold redshift. This is 
natural in the existing models of modified gravity that 
aim to explain the late-time acceleration, where depar- 
tures from GR occur at around the present day horizon 
scale. Also, the success in explaining the BBN and CMB 
physics relies on GR being valid at high redshifts. 

One could assume a functional parametrization for 
/i and 77, either motivated by a modified growth (MG) 
theory or by simplicity, and fit the parameters to the 
data [25, 32, 46]. We adopt a different approach, and 
pixelize p(a, k) and 77(0, k) on a grid in time and scale, 
fitting their values in each grid point to the data. We then 
solve the eigenvalue problem for the covariance of the pix- 
els (i.e. perform a PCA) to find their independent linear 
combinations that can be compared with their prediction 
in GR [30]. As we will elaborate later, the PCA method 
has several advantages, such as being model-independent 
and degeneracy-free, although it is much more compu- 
tationally expensive. In this work, we will utilize both 
functional fit and PCA strategies in order to search for 
any imprint of modified gravity. 

To begin with, we parametrize our Universe using: 



(uib, w c , 6 S , r, n s , A s , J\f, X) 



(19) 



where wt> = f2t>/i and uj c = Q c h are the physical baryon 
and cold dark matter densities relative to the critical den- 
sity respectively, & s is the ratio (multiplied by 100) of the 
sound horizon to the angular diameter distance at decou- 
pling, r denotes the optical depth to rc-ionization, and 
n s , A s arc the primordial power spectrum index and am- 
plitude, respectively. We also vary and marginalize over 



several nuisance parameters denoted by Af when perform- 
ing the likelihood analysis for weak lensing and SNe, as 
we will elaborate later. The modification of gravity is en- 
coded in X, and we consider two different kinds of MG 
parametrizations, Xj and Xu, as described in the fol- 
lowing subsections. Finally, we assume a flat Universe 
and an effective dark energy equation of state w = — 1 
throughout the expansion history. 



A. A single high redshift transition in /j, and ly. 

Xi = {7*0,770, Az} 

There is no physical reason to assume that depar- 
tures from GR ought to be scale-independent, and a 
PCA forecast analysis [30] actually showed that the scale- 
dependence of jji and 77 is better constrained than their 
average values or the time-dependence. Nevertheless, we 
shall first consider the case in which fi and 77 arc taken 
to be scale-independent and transit from their GR values 
to another constant value below a threshold redshift z s . 
Aside from simplicity, this will allow an insightful com- 
parison with the results of the scale-dependent analysis 
later. 

To model the time evolution of y, and 77 we use the hy- 
perbolic tangent function to describe the transition from 
unity to the constants /Jo and 770 as 

M( z ) = — tt^ (l + taim — \ — ~1 + Mo , 



■q(z) 



2 

1 - Vo 



1 + tanh 



Az 
z - Zt 
Az 



Vo 



(20) 



For a given z s , the above parametrization has three free 
parameters: /Joj'7o an d Az. In Figs. 1 and 2, we show 
the imprints of MG on the weak lensing aperture-mass 
dispersion (see details in Sec. II C) and CMB TT power 
spectra for different values of the MG parameters respec- 
tively. To be more physically transparent, we will also 
present our results in terms of the {/j, £} parametriza- 
tion (see Eq. (18)). 

We find that both WL and CMB obscrvables are more 
sensitive to the variation of £0 than to that of po for both 
the z s = 1 and the z s = 2 case. This is expected since WL 
and CMB (via the ISW effect) measure, respectively, the 
power spectra of ($ + \I r ) and their time derivative, and 
are primarily controlled by So- However, at late times on 
sub-horizon scales, po only affects ($ + ty) indirectly by 
altering the growth rate of A via 



A + HA - A^Gpa 2 p a A = , 



(21) 



where the over-dot represents differentiation with respect 
to the conformal time, and % = a/a is the conformal 
Hubble parameter. Therefore changing ($ + >!') by vary- 
ing /Jo is much less efficient than tuning the multiplier 
So- This can be seen in panels (Ci,C2), where we plot 
the relative difference of the evolution of A with respect 
to GR for different values of /jq- One can read from the 
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[arcmin] 9 [arcmin] redshift z 



FIG. 1: Imprints of modified gravity parametrized by X\ on the weak lensing aperture mass dispersion (panels Ai, A2), relative 

1/2 

difference in M a p with respect to GR (Bi, B2), and relative difference in A with respect to GR (Ci, C2). The model parameters 
are shown in the legend of panel (Ai). The shaded regions in panels {A\, A2, Bi, B2) are excluded from our analysis. The data 
with error bars over-plotted in panels (^1,^2) are taken from the CFHTLS survey. 
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FIG. 2: Imprints of modified gravity parametrized by X\ on 
the CMB TT power spectra for different threshold redshift 
z a and different transition width Az. Different models are 
distinguished by different line styles and colors, as shown in 
panel (Ai) of Fig. 1. The data points with error bars are taken 
from the WMAP 5-year survey. 



plot that, for example in the case z s = 2, if one fixes 

So = 1 and increases /j,q by 50%, A is enhanced by 15% 

at z = 0. From this, it follows from Eq. (18) that ($ + 

1/2 

is also enhanced by 15% at z = 0. Finally M a p should 
be enhanced by roughly the same amount according to 
Eq. (8), and this is what we see in panel (B2) (navy 



dashed-dotted line). On the other hand, if /Iq is fixed to 

unity and So is enhanced by 50%, then Eq. (18) clearly 

1/2 

shows that ($ + \f r ) and M a p should be enhanced by 50% 
as well, as shown in panel (B2) (magenta short-dashed 
line) . Therefore we can conclude that our obscrvablcs are 
more sensitive to E than /iq for both cases of the tran- 
sition redshift, z s = 1 and z s = 2. 

When varying the growth rate controller /zo, the ear- 
lier the redshift at which it is turned on, the more total 
change in growth and gravitational potential will be accu- 
mulated by the present day, as shown in panels (C\, C2). 
Thus our observables are more sensitive to the same 
amount of variation in the MG parameters in the case 
of z a = 2 compared to z 3 = 1 due to this "accumulation 
effect". 

In Fig. 2, we see that the CMB angular spectrum will 
strongly disfavor a sharp transition in S. This is obvious 
from Eq. (18) — a sharp transition in S triggers a sud- 
den change in (<& + *!>), making ($ + 4>) diverge around 
the transition region, and this is what we see in panels 
(Ai,Bi). Note that if the transition width Az is small 
enough, the ISW signal in the CMB spectrum converges 
because the details of the transition become irrelevant 
in the Az — > limit; then the ISW signal is determined 
solely by the change of ($ + \&) and the window function 
at the transition. We find that the choice Az = 0.05 is 
narrow enough to approximate well to this Az — > limit. 
We also show the results for a milder transition, namely, 
Az = 0.5 in panels (A 2 , B 2 ) of Fig. 2. We see that the (j,q 
curves arc largely unchanged, implying that the effect of 
varying (1q is not sensitive to Az. However, the 'bumps' 
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on the So curves become less pronounced for smoother 
transitions, as expected. In the following analysis, we will 
show results corresponding to the parametrization X\ for 
both fixed Az = 0.05, and a floating Az, which will 
then be treated as a nuisance parameter and marginal- 
ized over. 



B. Pixellation + PCA: X u = {^, (i = 1..4)} 

Even though the parametrization X\ has the advan- 
tage of simplicity, it is not theoretically well-motivated. 
Models of modified gravity commonly introduce a scale 
into the theory, and correspondingly \x and r\ always have 
some scale dependence. Moreover, this parametrization is 
not phenomenologically efficient to capture a deviation 
from GR. As shown in [30], the growth observables are 
much more sensitive to the scale dependence of the two 
functions, than to their time dependence. 

In order to be more general, one can pixelize \x and 
r\ in the (k,z) plane and treat their values on each grid 
point as free parameters. These parameters are in general 
correlated with each other, and this blurs the interpreta- 
tion if one attempts to constrain them directly. Instead, 
one can construct new variables that are uncorrelated 
linear combinations of the original parameters and use 
them to test GR. This can be achieved by diagonalizing 
the covariance matrix of the original pixels and using the 
de-correlation matrix to map the original pixels onto the 
uncorrelated variables. Such de-correlation, or PCA, has 
been used to study constraints on the evolution of the 
dark energy equation of state w(z) [47-49]. Here we em- 
ploy a two-dimensional PCA in the (fc, z) plane since /i 
and rj are functions of time and scale. 

The model-dependence disappears in the limit of a very 
fine tessellation of the two functions into pixels. In real- 
ity, computing costs limit the number of pixels one can 
afford to fit. To determine the optimal number of pixels, 
we performed a Fisher Matrix PCA forecast, analogous 
to the one in [30], finding that in order to capture the 
(k, z) dependence of the best constrained combined eigen- 
modes one needs at least a 2 x 2 pixels for \x and £ in the 
range of k £ [0, 0.2], z € [0, 2], as illustrated in panel A of 
Fig. 7. Note that here we choose to work with the {/i, £} 
parametrization because the ISW and WL constrain £ 
more directly than rj. We have checked that our results 
do not change if we pixelize {/i, 77} instead. To properly 
set the transition width between the neighboring pixels, 
we start from a wide transition width (Az = 0.5), and 
reduce it until the final results converge. As in the case 
of parametrization X\, we found convergent results when 
Az < 0.05, therefore we chose Az = 0.05 for the transi- 
tion width. 

Thus, in model Xu, we start by fitting 8 pixels, 
(i — 1..4)}, along with the non-MG parameters, 
to obtain the covariance matrix of all parameters. We 
then diagonalize the 8x8 block of the covariance matrix, 



C(>,£) corresponding to fi and S: 

Co,,s) = Wk- X W T . (22) 

The rows of the de-correlation matrix W are the princi- 
pal components [50], or eigenmodes, while the diagonal 
elements of A, i.e. the eigenvalues, are the inverses of the 
variances on the uncorrelated linear combinations of the 
original pixels. Namely, we use W to rotate the original 
parameters, denoted by the vector p, into new uncorre- 
lated parameters q defined as 

'/< 1 ' »<</'.- »'<•-• ( 23 ) 

3 3 

In GR one has q = 0, since p = 1, therefore we can test 
GR by performing a null test on q. By construction, the 
eigenvectors are orthogonal and the g's have uncorrelated 
errors given by the inverses of the eigenvalues. 

IV. RESULTS 
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FIG. 3: 68% and 95% C.L. contour plots for {/ko, ?7o} and 
{fj.0, So} for two different threshold redshifts: z s — 1 (up- 
per panels) and z 3 — 2 (lower panels). In both cases the 
transition width is fixed to Az = 0.05. From outside in, the 
shaded regions in cyan, yellow and blue illustrate the con- 
tours derived from the data of CFHTLS+CMB shift parame- 
ters, CFHTLS+WMAP5 and CFHTLS+WMAP5+ISW, re- 
spectively. For the contours shaded in the same color, the light 
and dark regions show the 68% and 95% C.L. contour respec- 
tively. In all cases, the SNe data are combined, and the priors 
of cosmic age, BBN and HST are applied. The star denotes 
the GR values. 



Given the set of cosmological parameters P in Eq. (19), 
we calculate the observables, which include the CMB 
temperature and polarization spectra, the CMB/galaxy 
cross-correlation (which we often refer to as ISW), the 
luminosity distance for SNe, and the WL aperture-mass 
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TABLE I: The mean values of Mo, ??o and Eo with 68% and 95% C.L. error bars for different models and for different data 
combinations. Note that '£ A , iZ','CMB', 'ISW and 'WL' are short-hands for WMAP5 shift parameters, full WMAP5 data, ISW 
data and CFHTLS data explained in the text, respectively. 





Z s = 1 


z s = 2 




Mo 




So 


Mo 


Vo 


So 


Az fixed 


CMB 
CMB+ISW 
WL+e A ,R 
WL+CMB 
WL+CMB+ISW 


-i n +0.11+0.40 
i - u -0.13-0.34 
n Q7 +0.09+0.37 
U - J ' -0.13-0.29 
n R o+0.65+1.36 
u -°°-0.45-0.57 

0.95 ± 0.24+° j4 
0.90 ± 0.21+ il 


i 1+0.51 + 1.0 
1 - ± -0.48-0.74 

1.2 ± 0.50+° f 7 

i 7 +3.2+6.2 
-■-■'-1.6-2.3 

1.2 ±0.60^0 ^ 

1.3±0.56lH! 


1.0 ±0.03 ±0.06 
1.1 ±0.028 ±0.055 
0.86 ±0.39 ±0.74 

1.0 ± 0.033lo;°g 
1.0 ±0.027 ±0.05 


i 1+0.16+0.37 
i ' 1 -0.17-0.31 
i n +0.15+0.35 

-■-■U-o.ie-o^s 

n t-o+0.92+1.17 
U.OO_ 38 _o .55 

0.87 ±0.1510:25 
0.84 ± 0.13io!2i 


n QC+0.11+0.62 
u - yu -0. 18-0.47 
fl QO+0.10+0.55 
u - yo -0. 17-0.45 
o i +3.0+5.7 
z - 1 -1.5-2.2 

1.4 ± 0.44+°;% 
1.4 ± 0.39l°; 8 69 


1.0 ±0.025 ±0.05 
1.0 ±0.024 ±0.05 
0.89 ±0.19 ±0.33 
1.0 ±0.03 ±0.05 
1.0 ±0.025 ±0.05 


Az float 


WL+CMB+ISW 


i -,+0.62+0.80 
L - L -0. 34-0. 45 


O.gsl?;^ 1 ^ 


n QzL +0.08+0.12 
"■■7^-0.14-0.32 


0.87 ± 0.12±g-_?| 


1.3 ± 0.35±g;|g 


1.0 ±0.03 ±0.06 
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FIG. 4: 68% and 95% C.L. contour plots for {mo, Vo} and 
{/j,q, So} for two different threshold redshifts: z s = 1 (upper 
panels) and z 3 — 2 (lower panels). In both cases the transition 
width is fixed to Az = 0.05. From outside in, the shaded 
regions in yellow and blue illustrate the contours derived from 
the data of WMAP5 and WMAP5±ISW, respectively. For the 
contours shaded in the same color, the light and dark regions 
show the 68% and 95% C.L. contour respectively. In all cases, 
the SNe data are combined, and the priors of cosmic age, BBN 
and HST are applied. The star denotes the GR values. 
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FIG. 5: 68% (dark shaded) and 95% C.L. (light shaded) con- 
tour plots for {fJ,o, So} and {fto, ?7o} for two different threshold 
redshifts: z s — 1 (upper panels) and z 3 — 2 (lower panels). All 
the constraints are from the combined data of ISW, WMAP5 
and CFHTLS. To obtain the front, green contours, the tran- 
sition width is fixed to Az = 0.05, while the blue contours 
on the back layers show the case of a floating Az, which is 
marginalized over. The dashed curves show the covered con- 
tour edges. The star illustrates the GR values. 



Parametrization X\ 



dispersion M ap , using MGCAMB. We then fit the available 
CMB, ISW, SNe and WL to observations using a modi- 
fied version of the Markov Chain Monte Carlo (MCMC) 
package CosmoMC 4 [51], based on Bayesian statistics. Our 
main results are summarized in Figs. 3-10 and Tables I 
-IV. 



http : //cosmologist . inf o/cosmomc/ 



Let us start with the single-transition parametrization 
X\ of the MG parameters. In Fig. 3, we show the 68% and 
95% C.L. contours of {^o, Vo} and {^o, S } for the cases 
of z s = 1 (upper panels) and z s = 2 (lower panels) for 
different data combinations. Here the transition width is 
fixed to Az = 0.05. We show contours derived from the 
CFHTLS data combined with the CMB shift parameters 
I A and R given in [38], CFHTLS plus full WMAP5 and 
CFHTLS±WMAP5±ISW. In Fig. 4, we also show con- 
tours derived from full WMAP5 and WMAP5 + ISW. 
All cases include the SNe data, and the cosmic age, BBN 
and HST priors. 
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\ \ 

FIG. 6: 1-D posterior distributions of Jio,rjo, and Eo derived 
from the joint analysis of ISW, WMAP5 and CFHTLS data. 
The black solid lines show the cases of sharp transition, i.e., 
Az = 0.05, while the red dashed lines illustrated the cases 
where Az varies and is marginalized over. The upper and 
lower panels are for the z a = 1 and z a = 2 cases, respectively. 
The vertical dashed lines illustrate ji = rjo = So = 1 to guide 
eyes. 



One thing that can be noticed from Figs. 3 and 4 is 
that the z s = 2 models, shown in the lower panels, are in 
general better constrained than the z s = 1 ones, shown in 
the upper panels. This is due to the 'accumulation effect' 
explained in Sec. Ill A. 

In Fig. 3, the largest cyan contours show that the 
CFHTLS WL data combined with CMB shift parame- 
ters are able to constrain Eo at ~ 20% level, but only 
weakly constrain jiq and rjo, since WL observables are 
directly sensitive to the variation in So as explained in 
Sec. Ill A. For the same reason, one sees little degeneracy 
in the {So, jio} plane. However, in the {no, 770} plane, the 
contours show a banana shape, which indicates that jio 
strongly anti-correlates with 770. This is understandable 
— one can increase no to enhance growth and thus but 
then 770 can be lowered to decrease <E>, leaving $ + un- 
changed. Also, as discussed in Sec. Ill A, the sensitivity 
of the observables to no is comparable to their sensitiv- 
ity to rjo , although the former is slightly larger than the 
latter. This is the reason why the degeneracy is visible. 
From Table I, we see that CFHTLS combined with CMB 
shift parameters favor slightly lower values of no and So 
than unity, but GR is still within the ltr level. 

The constraints become much tighter if one includes 
the full WMAP5 data, as shown in the yellow contours 
in Figs. 3 and 4. This is mainly because the WMAP5 
ISW-ISW ACFs strongly penalize abrupt changes in So, 
as we show in Fig. 2. The constraints get even tighter 
when the ISW-gal CCFs data are added, as illustrated 
in the innermost blue contours. From Fig. 3 and Table I 
we see that for the case of z s = 1 GR is fully consistent 
with the combined data. For the z s = 2 case, GR is also 
consistent, but is on the la edge, indicating that a model 
with a lower no would be slightly favored when WL data 
arc included. GR is always closer to the best-fit model 
when WL data are not used, as we can see in Fig. 4. 



If the transition width Az is allowed to vary in the 
range of [0.05, 0.5] one could expect a dilution of the 
constraints. We observe the result in Figs. 5 and 6, which 
show the contours and 1-D posterior distributions for MG 
parameters for the cases of sharp transitions and float- 
ing transitions for all the data combined. As we can see, 
marginalizing over a floating Az degrades the constraints 
on /io, rjo and S by roughly 150%, 50% and 300%, respec- 
tively for the z s = 1 case, but there is little degradation 
for the z s = 2 case. Again, we see that GR is a perfect 
fit in the z a = 1 case, while a model with a lower no 
is slightly favored in the z s = 2 case when all data are 
combined. 

This can be understood as follows. As shown in panels 
(Ai, Bi) in Fig. 2, for the z s = 1 case, a sharp transition 
in S produces a huge bump on CMB TT spectrum at 
£ < 70, where the cosmic variance dominates the error 
budget. If the transition is mild, the bump structure be- 
comes less pronounced as shown in panel A2, thus there 
is less tension with the CMB data, which in turn loosens 
the constraints on the MG parameters. However, for the 
z s = 2 case, the bump appears at £ < 150 on the CMB 
spectrum, where WMAP5 has precise measurements. It 
is true that relaxing Az reduces the bump feature some- 
what as illustrated in panel B2; however the constraints 
on the MG parameters cannot be diluted to a large extent 
due to the high quality CMB data at 70 < I < 150. 

Note that the constraint on S from the ISW effect is 
much tighter than the current WL constraints, due to 
the sensitivity of the ISW to the gravitational transition. 
This means that the ISW data provide valuable informa- 
tion on the time evolution of modified gravity parame- 
ters. 



B. Parametrization Xu 

We present results for the second parametrization X\\ 
with and without the inclusion of WL data, as summa- 
rized in Figs. 7, 10 and Tables II, IV respectively. Starting 
from the analysis based on the full data set, in panel (A) 
of Fig. 7, we show our S and ji pixcllation, and in panel 
(B), we show the 1-D posterior distributions of the eight 
S and jjl pixels. As we found for the parametrization Xi, 
the S pixels are in general better constrained than the ji 
pixels. We find that the constraints on all the pixels are 
consistent with the GR prediction except for that of S 3 , 
which is S 3 = 0.80 ± 0.12 ± 0.22 (mean value with 68% 
and 95% C.L. errors). This means that S3 deviates from 
the GR prediction at an almost 2a level. However, the 
correlation among all of the eight pixels blurs the naive 
interpretation of the seemingly 2a signal. Thus, we follow 
the PCA prescription explained in Sec. Ill B and obtain 
the linearly uncorrelated parameters, q's, using Eq. (23). 

The three best constrained eigenmodes are shown in 
panel (C) of Fig. 7 and are fairly well-localized. One can 
clearly see that they primarily depend on the S pixels, 
which are expected to be better measured by the ISW 
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FIG. 7: Panel (A): the E,/i pixellation used in this work; Panels (B) and (D): the 1-D posterior distributions of p and q, 
which denote the original pixels of E and /i, and the uncorrelated linear combinations of the original pixels, respectively. The 
likelihood distributions are normalized so that the area of each distribution is unity; Panel (C): the first three eigenfunctions, 
i.e. values of W relating p and q via Eq. (23). These are derived from the joint analysis of ISW, WMAP5 and CFHTLS data. 



TABLE II: Mean values, and 68% and 95% C.L. constraints 
of the original pixels (left panel) and the uncorrelated linear 
combinations of the pixels (right panel). All data sets are used. 
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and WL. In particular, the eigcnmodc corresponding to 
qi (i < 3) received the largest contribution from E, (i < 
3). From the eigenmodes we can deduce the following 
relations between the q's and the original pixels: 



Qi 

<12 



-1 

-1 
-1 



0.85E! -0.52S 2 



0.85 

0.52E! 



0.52 
0.85E 2 



= 0.0 ±0.02 ±0.04 



0.52 ±0.85 
E 3 = -0.17 ±0.06+°;^. 



0.0 ±0.05 ±0.10 

(24) 



is 



Namely, Ei is strongly degenerate with E 2 , while E3 
largely independent of Ei or E 2 ; 

One can understand this by realizing that the current 
ISW data can put stronger constraints on the MG pa- 
rameters than the current WL data, due to the sensitiv- 
ity of the ISW to any modification of growth at z < 2. 



FIG. 8: Best-fit aperture mass power spectra (M ap ) for differ- 
ent MG parametrizations shown in different colors and line 
styles. Black solid: parametrization An with all the pixels 
varying; Red dashed: parametrization An with only E3 vary- 
ing; Blue dashed-dotted: GR; Navy dashed-dotted-dotted: 
GR with a fixed fi m = 0.3. The data points with error bars 
show the CFHTLS data; the shaded region is excluded from 
our analysis. 



Future WL surveys (including the upcoming final results 
from the full CFHTLS survey) will feature higher signal- 
to-noise and will provide tighter limits on MG parame- 
ters [30]; however, at present, the ISW component domi- 
nates in our data combination, making the linear combi- 
nations of the E pixels on large scales (k < O.lh/Mpc), 
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TABLE III: The relative improvement on x 2 with respect to the ACDM model for different models. The A\ 2 is shown for 
different data separately. The mean values with 68% and 95% C.L. error bars of constraints on fi m and as are also shown. 
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FIG. 9: The 68 and 95% C.L. contour plots of {fi m ,£ 3 } (left 
panel) for parametrization Xu, and {f2 m ,£o} (middle and 
right panels) for parametrization Xj. See text for details. 



Si and E 2 , best measured. 

The constraints on the q's are summarized in the right 
panel of Table II, and the 1-D posterior distributions are 
shown in panel (D) of Fig. 7. We find that all the g's are 
consistent with zero as predicted by GR, except for q^, 
which deviates from the GR prediction at more than 95% 
confidence level. This means that our measurement of E 
and/or fi deviate from unity at the level of at least 95% 
C.L. at some point in {k, z). 

To spot the source of this signal, we fix all the pixels to 
unity except for E3 , which we allow to vary. Interestingly, 
we find that the goodness-of-fit of this one-pixel model 
is almost identical to that of the 8-pixel model. To be 
explicit, we list the \ 2 f° r each dataset separately for both 
models, and also show the constraints on fi m and as for 
these models in Table III. For comparison, we also do the 
same for the models parametrized using X\. Comparing 
to ACDM, we find that allowing E3 to vary can reduce 
the WL x 2 by 2.2, while also reducing the SNe x 2 by 
roughly 12. This E 3 CDM model has a best-fit Q m of 0.3, 
which is much larger than 0.24 for ACDM. The allowance 
for a high f2 m is the reason for the significant SNe data 
preference for this one-pixel model, since n m = 0.3 is 
the best-fit value for the SNe sample we use [4]. Note 
that in ACDM, f2 m = 0.3 is strongly disfavored by WL 
data, as we show in Fig. 8 (navy dashed-dottcd-dottcd 
line) , since there increasing il m shifts the best-fitted M ap 



(blue dashed-dotted) on all scales, which is in serious 
disagreement with the data, especially on scales 8 < 60 
arc min. 

The cause of the apparent 2a hint of departure from 
GR can be easily identified. There is a clear "bump" fea- 
ture in the CFHTLS data (e.g. Fig. 8) at 8 ~ 120 arc 
min, which can be attributed to a systematic effect [34] : 
according to the CFHTLS team, this is a known issue, 
due to residual ficld-to-ficld variations in shear estima- 
tion on the scale of the camera field-of-view. As an in- 
formative exercise, we study how we could improve the 
fit assuming a cosmological source for the feature. One 
could shift the curve at 8 > 60 arc min to follow the 
"bump" more closely. Such a scale-dependent tweak of 
Map can not be realized by tuning the MG parameters 
in the parametrization X\. However, in the parametriza- 
tion X\\, one can achieve this by firstly increasing Q m , 
then lowering the growth rate on small scales (8 < 60 
arc min), which can be effectively done by lowering E3. 
The resultant fit is shown in Fig. 8 as a red-dashed line, 
which is almost identical to the best-fit 8-pixel model. 

To see this point more clearly, in Fig. 9 we show the 
contour plots between fi m and E3, and fl m and Eo of 
X\ for the combined data. As we can see, in the one- 
pixel (E3) model can one obtain a high il m as favored by 
the SNe data, while in parametrization X\, Eo is tightly 
constrained, so that a high f2 m is definitely not allowed. 
Notice that the CMB and ISW data disfavor a high f2 m , 
as we see in Table III, but the preference from the WL 
and SNe data outweighs this penalty, making a lower E 3 
and higher f2 m strongly favored by the combined data. 
Also, the one-pixel model would be strongly favored by 
the data from the model-selection point of view, since one 
can reduce the total x 2 by 8.6 by introducing one more 
parameter (E 3 ) over ACDM, even though this model was 
constructed a posteriori. 

As stated above, it is likely that the "bump" , which is 
responsible for this 2a deviation, is due to a systematic 
effect [34], which according to the CFHTLS team is due 
to residual field-to-field variations in shear estimation on 
the scale of the camera field-of-view; this explains the 
scale of the bump. On these grounds, we stress again that 
it is premature to make any statements about the validity 
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of ACDM based on this feature, even though technically 
we cannot rule out the new physics at this point. 

In order to be conservative about this issue, we also 
perform the X\\ analysis without including the WL data 
from CFHTLS. The results are summarized in Fig. 10 
and Table IV. The relation between the uncorrclated pa- 
rameters g's and the original pixels, and the 68% and 
95% C.L. constraints on the g's are now given by 
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As expected, the first two best constrained modes are al- 
most unchanged even if we remove the CFHTLS data, 
confirming that these modes are mostly constrained by 
the ISW effect. On the other hand, the bound on S3 be- 
comes very weak, demonstrating that the CFHTLS data 
are responsible for the constraint on this parameter. In 
the current case we find that all the eight uncorrelated 
parameters are consistent with GR with 95 % C.L. 
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FIG. 10: Panels (A) and (B): The 1-D posterior distributions 
of p and q, which denote the original pixels of E and n, and 
the uncorrelated linear combinations of the original pixels, 
respectively. The likelihood distributions are normalized so 
that the area of each distribution is unity. These are derived 
from the joint analysis of ISW and WMAP5 data, without 
the CFHTLS WL data.. 



V. CONCLUSION AND DISCUSSION 

We have tested GR with current cosmological data, 
using a framework in which the departures from GR are 
encoded as modifications of the anisotropy and Poisson 
equations; these equations specify, respectively, how the 
metric perturbations relate to each other, and how they 
are sourced by perturbations in the energy-momentum 
tensor of matter. The modifications were parametrized 
with two functions {77, /it} (or alternatively {E,/it}) that 
reduce to unity in GR. We have then explored the con- 
straints on these functions in two ways. First, we have 



TABLE IV: Mean values, and 68% and 95% C.L. constraints 
of the original pixels (left panel) and the uncorrelated linear 
combinations of the pixels (right panel). Results without the 
CFHTLS WL data. 
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allowed them to evolve from unity at high rcdshifts to 
a different value today in a scale-independent way. Sec- 
ond, we have pixelized them in both scale and redshift and 
performed a Principal Component Analysis (PCA), fol- 
lowing the ideas of [26, 30] — a first general study of this 
kind. Specifically, we have used a 2 x 2 pixellation for each 
function, thus having 8 modified gravity parameters. In 
order to remove the covariance between the bins, and to 
analyze which modes are best constrained, we have then 
performed a 2D PCA of the results, obtaining constraints 
on the 8 derived de-correlated parameters. 

We have used currently available data constraining 
both the background expansion history and the evolu- 
tion of scalar perturbations in the Universe. In particu- 
lar, we have used a combined measurement of the ISW 
effect through correlation of galaxies with CMB, the lat- 
est available supernovae Type la data including those 
from the SDSS, the CMB temperature and polarization 
spectra from WMAP5, and weak lensing data from the 
CFHTLS shear catalog. We have kept the analysis con- 
servative by excluding small-scale data in the strongly 
non-linear regime, and we have checked and excluded 
possible tensions between the data sets by analyzing 
them individually before combining them. 

Throughout the paper, we have assumed a flat ACDM 
background and tried to constrain deviations from GR 
in the evolution of matter and metric perturbations. In 
the simplest case, where the MG functions {/z, £} were 
allowed a single transition in redshift, we have found no 
evidence for a departure from GR, in agreement with 
other works. We find that the ISW effect, probed through 
the CMB auto-correlation and the cross-correlation with 
galaxy maps, currently gives the strongest constraint on 
£ because it is sensitive to the change of the lensing po- 
tential, $ ± ^, at the transition. 

In the pixellated case, we have found that one of the 
PCA eigenmodes shows a 2a deviation from GR. How- 
ever, this anomalous mode is due to the "bump" feature 
in the CFHTLS lensing data, which is most likely due to 
a systematic effect [34], combined with a preference for 
higher Q m by the SNe data. A separate analysis which 
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does not include WL data shows indeed good agreement 
with GR. A better understanding of systematic effects in 
both WL and SNe data sets needs to be achieved before 
any such discrepancy is viewed as a deviation from GR. 

Even though this is most likely due to a known system- 
atic effect, we emphasize that we would not have found 
this deviation if E were taken to be scale-independent. 
In such case, the change in E would be significantly con- 
strained by the ISW effect. The PCA analysis using two 
bins in k for E could successfully isolate the strong con- 
straint from the ISW effect and pick up a feature in 
WL. This demonstrates that the same data can lead to a 
higher level of detection of deviations from an expected 
model if more flexibility is allowed in the parametriza- 
tion. Thus, when fitting /i(k, z) and r](k, z) to data, it 
is important that their parametrization allows for some 
scale-dependence. Otherwise, one might risk missing a 
systematic effect or a real departure from GR, and thus 
would not be exploiting the true discovery potential of 
the data. 

Finally, we comment on other recent studies that re- 
ported constraints on deviations from GR using current 
cosmological observations. In [33], the COSMOS weak 
lensing tomography data [9] were used together with SNe, 
CMB, BAO and the ISW-galaxy cross correlation. [32] 
found that the constraints from CMB+SNe+CFHTLS 
without COSMOS were indistinguishable from those in- 
cluding COSMOS and they did not find any deviations 
from GR. They argued that the sky coverage of CFHTLS 
is more important than the redshift depth of COSMOS. 
Also it should be noted that weak lensing measurements 
in COSMOS are made on strongly non-linear scales and 
there is an ambiguity in modeling the non-linear power 
spectrum. 

In [32] , a similar set of data to that described here was 
used to constrain two functions that are combinations of 
(i and rj. The differences between our study and that of 



[32] include: (1) we used the ISW cross-correlation data; 
(2) we excluded small-scale modes in the CFHTLS data 
to avoid the strongly non-linear regime; (3) we simul- 
taneously constrained two functions fi and rj while [32] 
varied only one of the parameters when they use 3 bins 
in z; (4) scale dependence was allowed in our paper; and 
(5) we used the Fisher matrix based PCA approach to 
make a decision on how many pixels to use. Their results 
are qualitatively consistent with the result of our first 
parametrization. 

Further improvements of this technique will be pos- 
sible with a new generation of LSS data (e.g. DES, 
Pan-STARRS, LSST, Euclid), which will dramatically 
increase the number of modes with sufficient signal-to- 
noise. Finally, peculiar velocity data will provide an ad- 
ditional valuable probe for our approach, since they can 
constrain [i directly [52] , thus breaking the H-fj, degener- 
acy. 
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